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Abstract 

Renormalization group method is applied to the study of vibrational en- 
ergy transfer in protein molecule. An effective Lagrangian and associated 
equations of motion to describe the resonant energy transfer are analyzed in 
terms of the first-order perturbative renormalization group theory that has 
been developed as a unified tool for global asymptotic analysis. After the 
elimination of singular terms associated with the Fermi resonance, ampli- 
tude equations to describe the slow dynamics of vibrational energy transfer 
are derived, which recover the result obtained by a technique developed in 
nonlinear optics [S.J. Lade, Y.S. Kivshar, Phys. Lett. A 372 (2008) 1077]. 
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1. Introduction 

The study of vibrational motions in protein system has attracted much 
interest in the context of their relations to biological functions, paying at- 
tention to the intrinsic flexibility of large biomolecules. A lot of theoretical 
investigations on the basis of the normal mode or principal component anal- 
ysis [l-5] have been carried out in combination with molecular mechanics or 
dynamics simulations for protein systems, thus providing useful information 
at molecular level. Vibrational energy transfers between the normal modes 
are then expected when the nonlinear, anharmonic couplings are taken into 



accout 



BEi. 



The dynamical relaxation behavior after the excitation of vi- 
brational modes is thus central for a deep understanding of the relationships 
among structure, dynamics and function in protein systems. 



In this context, Moritsugu et al. 



9( studied the intramolecular transfer 



of vibrational energy in myoglobin by means of molecular dynamics simula- 
tion. They found that the vibrational energy was transferred from a given 
normal mode to very few other modes that were selected by the frequency 



relation of the Fermi resonance [10| and the magnitude of the third-order 
mode coupling term. It was also observed that the magnitude of the coupling 
coefficients was mainly determined by the degree of geometrical overlap of 
the corresponding normal modes. Further, performing a molecular dynamics 
simulation with just these few modes, they reproduced the results for all- 
atom simulation well, thus confirming that only these modes were important 
for the dynamics of vibrational energy transfer. 

As is well known, the description by naive perturbation theory breaks 
down for this kind of resonant energy transfer problem due to the emer- 
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gence of singular terms associated with the Fermi resonance 10j. Invoked 
by a theoretical technique for describing the resonant parametric interaction 



between different harmonics in nonlinear optics, Lade and Kivshar [11| pro- 
posed a simple theory to overcome this difficulty, which provided coupled 
mode equations for slowly varying envelopes of the vibrational modes. They 
thus found that their equations could describe the dynamics of the resonant 
energy transfer between the vibrational modes observed in the simulations by 
Moritsugu et al. js, 9] with good accuracy, and gave a wealth of suggestions 
about the underlying mechanism of nonlinear dynamics of vibrational energy 
transfer [12]. 

In this work, the renormalization group (RG) approach which has been 



developed as a unified tool for global asymptotic analysis |13Ml5| is employed 
in order to describe the long-time dynamics of the resonant vibrational en- 
ergy transfer in protein system. Starting with the same model Lagrangian 
as studied by Moritsugu et al. 9} and Lade et al. [llj], we apply the RG 
technique at the first-order perturbation level and derive the corresponding 
amplitude equations through the elimination of singular terms associated 
with the Fermi resonance. We then find that the derived equations are es- 
sentially of the same form as those obtained by Lade and Kivshar 11| at 



the lowest-order level. In the following section, the theoretical formulation 
of the present RG approach is addressed. Some discussions and conclusion 
are given in Sec. 3. 



2. Theory 



We start with a model Lagrangian given by Moritsugu et al. 
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4 1 4 1 

L = 2^ ~ ~ aqiq2q3 ~ W 

i=i i=i 

where g.,- refers to the coordinate of the jth normal mode with the frequency 
Uj and the dot denotes the time derivative. Two types of third-order coupling 
among the modes are considered in the model, which are characterized by 
small coupling parameters a and /3. Equations of motion for qj then read 

q\ = -u\qi - aq 2 q 3 - 2/3gig 4 , 



q 2 = -u\q 2 - aq x qz, 



q 3 = -u1q 3 - aq x q 2 , 



q\ = -ulq 4 - (3qf. (2) 
These equations of motion were found to appropriately reproduce the tempo- 



er observed in the all-atom molec- 



ral variations of vibrational energy trans : 
ular dynamics simulation for myoglobin 

The solution to qj is given by a perturbation theory up to the first order 
of a and f3 as 

m (t) - AiP iuJlt - aA ^ A 3 i(u>2+w 3 )t _ aA*A 3 M^-^t 
qi{i — i\\e 2 e 2 e 

2/3AtA 4 j(uj l+W4 )t _ 2(3AIA 4 jj^-^t , 



Uj\ — {bj\ + UJ4) 2 uj\ — (w 4 — UiY 
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q 2 (t) = A 2 e 1 ^ - 



aA x A 3 
u 2 - (wi + u 3 ) 2 



aAiAl 
lo\ - (uji - u 3 ) 2 



3 i(u>i-w 3 )t 



+ C.C. 



q 3 (t) = A 3 e^ 



uA x A 2 



aAiAl 



u 2 - (ui - u 2 y 



J(u 1 -U)2)t 



+ C.C, 



q 4 (t) = A^ - J^—^ - ^1 4 cc. 



uj{ - Auj{ 



^'7 



(3) 



where we have introduced amplitude parameters, Aj, and cc. (or the asterisk) 
means the complex conjugate. 

Here, we consider a case of vibrational energy transfer in which the Fermi 
resonance takes place. According to Moritsugu et al. jsl, 0], it is assumed 
that two frequency parameters, Q a = ui—u 2 — u 3 and ftp = 2uj\ —U4, satisfy 
such a condition as \Q a \, \£lp\ *C u)j. Retaining only the terms of the first 
order of a or /3 that become large under the Fermi resonance condition on 
the right-hand side of Eq. (3), we find 



aA 2 A 3 r _ inat _ 2(5 A\A^ 



-iQat 



ft a (uJi + UJ 2 + UJ 3 ) 



e Wlt + cc 



92 (t) 



A 2 + 



aA 1 An 



+U 2 - UJ 3 ) 



e iLU2t + c.c, 



q 3 (t) 



3 fl a {ui -u 2 + UJ 3 ) 



e wzt + cc, 



?4(*) 



.4, 



0AI 



Vtp{2uji + to 4) 



e iuJ4t + cc. 



(4) 
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According to the RG technique 13Ml5l|. we express the amplitude param- 
eters Aj in Eq. (4) up to the first order of a or (3 as 



A 1 = A 1 (r)[l + aZ la (r)+(3Z 1 ^r)], 
A 2 = A 2 (r) [1 + aZ 2a (r)} , 
A 3 = A 3 (r) [l + aZ 3a (r)\, 

A 4 = A 4 (r)[l + /3Z^(r)}, (5) 

where we have introduced a time variable r and renormalization constants 
Zj a , Zjp. In order to eliminate the singular terms in Eq. (4) up to the lowest 
order, the renormalization constants are required to satisfy the following 
equations: 



Ai{T)Z la (T) - 



A 2 (r)Z 2a (r) 



A 3 (r)Z 3a (r) 



Q a (oji + u 2 + oj 3 ) 
2i*(r)i 4 (r)e-^ r 

VlfjUJ 4 

ii(r)i*(r)e^« T 
Q (wi + u 2 -u 3 )' 

ii(r)i^(r)e^ ar 
Q a (ui -u 2 + u 3 ) ' 
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i 4 (r)Z 4/3 (r) 



AArfe 



(6) 



Since the amplitude parameters A, should not depend on the introduced 
time variable r, we call for identities as 



IF- ' 



(7) 



We then obtain the RG equations 
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15| up to the first order of a or (3 as 



dii(r) 
dr 



+ 



dr 



dr 



A 1 (r)=0, 



dA 2 (r) dZ^jr) - 

— h a — A 2 [t) = 0, 

dr dr 



dA 3 (r) dZsJj) - 



(8) 



dr dr 

from Eq. (5), considering dAj{r)/dr ~ 0(a) or 0(j3). We also note the 
lowest-order relations as 



AAr 



dr 



iA 2 (T)A 3 (r)e' in - T 

UOi + U0 2 + 



dZ 1/3 (r) _2ii*(r)i 4 (r)e-^ r 
dr 



Afr 



dZ 2a (r) _ lA^A^ry^ 
dr U\ + uj 2 — u 3 
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~ dZ^r) zi 1 (r)i*(r)e^- 
A s{V — = ■ , 

GST UJi — OO2 + UJ3 

Ai(T) ^ij) = JMrfe^ 
K J dr 2^1 + ^4 

from Eq. (6). Combining Eqs. (8) and (9), we finally obtain differential 
equations, 



dAi(r) iaA 2 (r)A 3 (T)e^ naT 2i/3A*(r)A 4 (r)e-^ r 



dA 2 (r) _ iai 1 (r)i*(r)e i ^ T 



dr loi + L02 — 

dA 3 (r) _ 2ai 1 (r)i^(r)e^ T 
dr u\ — (jJ 2 + 0J3 

dA 4 {r) if3A 1 (T) 2 e i ^ T 



(10) 



dr 2u>\ + U4 

to determine the renormalized amplitudes Aj(r). 

We now equate the time variable r with time t in Eq. (10), and solving 
the differential equations with the initial condition of qj(0) = Aj(0) + A*(0), 
we can obtain Aj(t). Thus, we find a solution to qj as 



q 3 {t) = Ajfte*** + i*(t)e-^'*. (11) 

v of each mo 
tion BUlill. 



The energy of each mode can then be evaluated in the harmonic approxima- 
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3. Discussion and Conclusion 

In the preceding section we have derived an amplitude equation, Eq. (10), 
on the basis of the RG approach. This equation is free from the singularities 
associated with the Fermi resonance, which have been observed in the naive 
perturbation th eory , and describes a nonlinear, slow dynamics of vibrational 
energy transfer [12| in terms of the renormalized amplitudes Aj. The present 
renormalization approach has thus succeeded in an effective separation of the 
two time domains associated with the fast dynamics governed by Uj and the 
slow ( "envelope" ) dynamics for vibrational energy transfer in protein. 



Lade and Kivshar |ljj analyzed the same model as in the present study, 
and derived an analogous set of amplitude equations by integrating out the 
fast vibrational motions in an intuitive way developed in the field of nonlinear 
optics. By solving the resultant equations numerically, they found that the 
long-time behavior of the vibrational energy transfer in the model Eq. (1) can 
be described very well, thus demonstrating the excellence of their method. 
Taking into account the condition that \Q a \ and \Qp\ <C Uj, their derived 
equations are essentially identical to the present Eq. (10). That is, their 
result can be recovered in terms of the first-order perturbation theory in 
the RG approach, and can be systematically improved in the more general 
framework of the present theory. 

In order to explain the dynamical behavior of vibrational energy trans- 
fer observed in their all-atom molecular dynamics simulation for myoglobin, 
Moritsugu et al. 0,8 considered a model system given by Eq. (1). Perform- 
ing a molecular dynamics simulation for the model system with the param- 
eter set of wi = 5.385, u 2 = 2.362, u 3 = 3.024, u 4 = 10.779, a = -0.13 and 
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(5 = 0.12, they reproduced the temporal evolution of the resonant energy 



transfers among the normal modes quantitatively. Lade and Kivshar 
then solved their amplitude equations numerically with the same parameter 
set as that by Moritsugu et al., and found a good agreement between the two 
results for the time course of energy variations of the normal modes. The en- 
ergy transfer dynamics observed in their studies is essentially a renormalized 
one in the sense that a slow dynamics emerges in contrast to the normal- 
mode dynamics characterized by Uj. The present work reproduces the result 
by Lade and Kivshar at the lowest-order level, and provides a more general 
framework for describing the resonant vibrational energy transfer in protein 
on the basis of the RG theory. This line of theoretical development would 
thus give a key tool for performing a coarse-grained description of protein 
dynamics in a comprehensive way. 
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